Decoding the riddle of F(4260) and Z c (3900) 
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The observation of Z c (3900) by the BESIII collaboration in the invariant mass spectrum of J/i^-k 
in e + e~ — > J/ipTr + -K~ at the center of mass 4.260 GeV suggests the existence of a charged DD* +DD* 
molecular state with I(J F ) = which could be an isovector brother of the famous X(3872) and 

an analogue of Zf,(10610) claimed by the Belle Collaboration. We demonstrate that this observation 

£C) ' provides strong evidence that the mysterious Y(4260) is a _DDi(2420) + -D-Di (2420) molecular state. 

7-H . Especially, we show that the decay of this molecule naturally populates low momentum DD* pairs 

and leads unavoidably to a cusp at the DD* threshold. We discuss the signatures that distinguish 

£N| ' such a DD* cusp from the presence of a true resonance. 

5—( < 

C$ . PACS numbers: 14.40. Rt, 13.75.Lb, 13.20.Gd 

\^ t During the past years, the experimental observation of a large number of so-called X, Y, Z states has initiated 
. tremendous efforts to unravel their nature beyond the conventional quark model. Especially, the conformation of 
signals in charged channels would be a direct evidence for exotic states. For instance, the Belle Collaboration reported 
signals for Z(4430) in 'ijj'n ± , and Zi(4050) and ^(4250) in Xci^ m B meson decays However, an enhancement 
in the same mass range was interpreted as a reflection by BaBar 0] ■ The more recent experimental results for charged 
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i ■ bottomonium states -££,(10610) and Zt,(10650), located close to the BB* and B* B* thresholds, respectively, by the 



Oh 

CD 



Belle Collaboration Q seem to be the first strong evidence for QCD "exotics" in the heavy quark sector. In this 
context the recent report of an enhancement in the J/f/'7r ± invariant mass distribution around 3900 MeV, right at 
the DD* threshold, by the BESIII collaboration :JJ clearly reinforces the existence of such an unusual phenomenon. 
This state, called Z c (3900) below, might be an isovector partner of the well established 1 ++ isoscalar V(3872) 
; but with I G {J PC ) = l+(l+~) for the neutral state. 
I/-) i In this work we demonstrate that, if V(4260) is a DD\ + c.c. molecule (below we use DD\ as a short notation), 
' the appearance of an enhancement around 3900 MeV in the J/ipK invariant mass distribution can be shown to be 
natural. Here D x refers to the narrow axial vector Dx(2420) (r = 27 ± 3 MeV) with I(J P ) = |(1+) 0|. In this sense 
the observation of the charged Z c (3900) state by BESIII in y(4260) — > J/tpTTir provides a very strong evidence for the 



molecular nature of V(4260). We also discuss whether the observed enhancement can be interpreted purely as a cusp 
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or whether the inclusion of explicit poles in the DD* system is necessary. Before we go into details of our calculations 
we first briefly review the status of V(4260). 

The V(4260) was observed via its decays into J/ipTr + ir~, J/^7r°7r°, and J/ipK + K~ by various experiments. Its 
mass and total width are well determined as m = 4263lg MeV and Ty = 95 ± 14 MeV, respectively The new 
data from BESIII confirms the signal in V(4260) — > J/ipir + ir~ with much higher statistics. The most mysterious fact 
about V(4260) is not that its mass does not agree to what is predicted by the potential quark model. Instead, as a 
. . .' charmonium state with J PC = 1 , it is only "seen" as a bump in its two pion transitions to J/tp, but not in any 
open charm decay channel like DD, DD* + c.c, D*D* and D S D*, or other tens of measured channels. In fact, the 
cross section lineshapes of the e + e~ annihilations into Z)W meson pairs appear to have a dip at its peak mass 4.26 
GeV instead of a bump. 

In the vector sector with J PC = 1~~, one should recognize that the DDi(2A20) is the first open charm relative 
5-wave channel that is coupled to this quantum number, and the nominal threshold is only about 29 MeV above the 
location of the V(4260). The possibility that V(4260) may be a bound system of DD\ and D*Dq was investigated 
in Refs. 0-01 ■ It should be stressed, however, that broad components (the widths of Dq and Di(2430) are as large 
as 300 MeV) can not produce narrow resonances [10(. We thus do not consider the D*Dq or Z?*Di(2430) component 
here, but focus on the assumption that the V(4260) is a bound system of only £>Z?i(2420). 
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FIG. 1: The Feynman diagrams for Y(4260) — > DD\ + c.c. — > J/ipmr and h c Trn considered in this work. 



It should be emphasized that in the literature many other solutions were proposed for the y(4260) (see Ref. [11 1 
for a recent review). From a fit to the 7T7t spectrum of F(4260) — > J/tpirir from the BaBar Collaboration 12j, Dai 
et al. [l3j] concluded that the even lower threshold XcO 1 -^ would have the largest coupling to F(4260), while the DD\ 
coupling to F(4260) turned out to be negligible. However, the observation of the enhancement at the DD* threshold 
in the J/ilm invariant mass spectrum actually rules out such a scenario and suggests that the underlying dynamics 
should be more sensitive to the DD\ threshold. 

In this work we do not try to identify the mechanisms that lead to the formation of the y(4260) as a molecular 
state, but study the consequences of the assumption that this state is a DD\ bound system. We argue that the 
interpretation of Y (4260) as a relative S- wave DD\ system is able to accommodate nearly all the present observations 
for y(4260). Especially, its absence in various open charm decay channels mentioned above and the observation of 
Z c (3900) in y(4260) — > J/tfmir can be naturally understood. 

The heavy quark spin symmetry also implies the presence of a D*D\ and D* Di component in the wave function. 
However, they will be neglected here since their corresponding thresholds are almost 200 MeV above the mass of 
y(4260). Meanwhile, we note that in a more microscopic treatment there should be heavier spin partners of F(4260) 
that should contain the mentioned constituents prominently. We will briefly come back to this issue at the end of this 
paper. 

We now provide the details of our calculations based on the picture described above. The pertinent diagrams are 
shown in Fig. [T] The central goal of our study is to pin down the structure of y(4260) and identify the quantitative 
importance of the DD* cusp from the loop diagrams (b)-(d) in order to understand whether the BESIII spectra call 
for the additional inclusion of an explicit pole diagram as depicted in diagram (a). 

It should be stressed that the vicinity of the DD\ threshold to 4.26 GeV favors the formation of low momentum 
DD* pairs, which may lead to the formation of the Z c (3900) bound systems. The reason is that the DD\ intermediate 
state as well as the DD* intermediate state can be simultaneously close to their mass shells — all with relative S 
waves. This gives rise to a triangle singularity studied in a different context in Refs. [l4, 15 1. Such a two-cut condition 
strongly enhances the corresponding matrix elements. In this sense the DD\ intermediate system provides an ideal 
doorway state for a low momentum DD* system. There are a series of S-w&ve open charm thresholds with J PC = 1~~ 
around 4.26 GeV, i.e. L>£>i(2420), Z)L>i(2430), D*D , D s0 D*, D s D sl , and £)*L»i(2420). However, all of them, except 
for Z)Z?i(2420), are either far away from the observed physical F(4260), or too broad to make a bound state [5j. Thus, 
we do not take them into account explicitly here. 

Note that in order to formulate the problem as an effective field theory the power counting of Refs. [l6|, [TtJ needs 
to be adapted to the present situation. We leave this to be reported in a subsequent work and focus here on a more 
phenomenological investigation. 

A complete calculation would also need the inclusion of the 7T7T final state interaction. This should be done in a 
way consistent with the Watson theorem. Thus, an inclusion of final state interactions in all diagrams is needed. In 
addition, since the pion pairs are in the isoscalar channel, also KK intermediate states need to be considered. This 
is beyond the scope of this paper and will be studied in a followed-up work. However, as a result of the omission 
of the 7T7r final state interaction we do not expect to be able to properly describe the 7T7t invariant mass distribution 
observed by BESIII. 
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FIG. 2: The invariant mass spectra for (a) J/ipn and (b) nn in Y(4260) — > J/ipTT-w. The solid, dashed, and dot-dot-dashed lines 
corresponds to the results of the full calculation, box diagrams, and triangle diagram with the Z c (3900) pole, respectively. 



Since Z c (3900) has the same quantum numbers as Z b most interactions needed for this work can be taken from 



Ref. [17|] . Here we only present the interactions between a P-wave charmed meson and other fields since they play 
a crucial role in this work. The heavy quark spin symmetry allows heavy mesons to form spin doublet super fields 
distinguished by their light degrees of freedom s\ — s q + I with s a the light quark spin and I the orbital angular 
momentum. For I — 1 it can be classified into two super fields 191. The interaction terms relevant for this work 
read 



C Y = iy(DtY i D i ? a -D? a Y i Di)+H.c, (1) 
6D^(d'd^ ab )Dl b + D*J\d^d^ ab )D\ b ] + H.c, (2) 



C Dl = i^[&D\ a {d l d^ ab )D*^ -D\ a {d^d^ ab )D^ 



where D (D') and D (D^) contain the annihilation (creation) operators for cq and cq fields, respectively. The analogous 
conventions are applied to D* and D±. The interactions contain five coupling constants in total. Since we focus on 
the shape of the invariant mass distributions only, all couplings that are common to all diagrams are not relevant. In 
this sense the only free parameters that influence the invariant mass distributions are the mass of Z c (3900) and the 
coupling constant z for Z C DD*. 

In analogy to what is known about the Z b states we assume that the width of the Z c (3900) is saturated by the DD* 
intermediate state. Then the width of Z c that enters the Flatte distribution is controlled by the same parameter z 
introduced above. 

For simplicity in this exploratory study we treat V(4260), £>i(2420), D and D* as stable states. Their masses are 
taken from the Particle Data Group @. For the Z c (3900) propagator we use a Flatte distribution, as mentioned 
above. By assuming that V(4260) is dominated by the Z)Di(2420) molecule component, we can estimate the YDD\ 
coupling by Weinberg's compositeness theorem 0, [2l|, i.e. y 2 /47r = 4(m£> + m^J 5 / 2 ^JlbE jmumu 1 — 17 GeV 2 . 
This predicts the dominant decay of V(4260) — > DD*tt + c.c. via the intermediate state l)Z?i(2420) + c.c. of which 
the partial width is larger than 40 MeV. This value is consistent with the total width measured for V(4260). In 
addition, the DD\ component should also distort the invariant mass distributions for F(4260) — > DD*ir + c.c. A more 
quantitative result will be reported in a subsequent work. 

The numerical results for the J/ipir and irir invariant mass spectra from F(4260) — t J/ipirir are shown in Fig. [21 
The results from the box diagrams (Fig. Q] (b)-(d)) are shown as the dashed lines, while the result for the Z c (3900) 
pole diagrams, Fig. Q] (a), is shown as the dot-dot-dashed lines. The full result, from the interference of these two 
contributions, is shown as the solid lines. One can see that even without the explicit inclusion of Z c (3900), two 
structures appear at the same masses in the J /^nr spectrum as in the BESIII data. Note that by charge conjugation 
invariance the J/tpir~ spectrum is identical with that for J/ipir + . The cusp at Mrj + Mjj* ~ 3.876 GeV in Fig. dja) 
marks the DD* threshold. The lower bump between 3.4 and 3.6 GeV in the J/i/>7r + invariant mass distribution comes 
from the interference between the conjugate diagrams where either a ir~ or a ir + is emitted first. Therefore, the lower 
bump is simply a reflection of the narrow structure at 3.9 GeV. 

Although one also observes that the agreement of the dashed curve in Fig. [2] with the data is only qualitative, 
the appearance of a visible cusp supports the claim made in the introduction that the loop kinematics favors the 
production of low momentum DD* pairs. However, the effect of the cusp alone is too small and too narrow compared 
to the data. This implies that the experimental observation of a very pronounced, broader peak at the DD* threshold 
suggests the need for an explicit resonance pole on top of the cusp. In this sense, the BESIII results seem to have 
confirmed the existence of the isovector partner of V(3872). 
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FIG. 3: The invariant mass spectra for (a) h c n and (b) -kit in Y(4260) — > h c nn. The notations are the same as in Fig. [2] 



The solid lines in Fig. [2] show the result of the sum of all diagrams in Fig. [TJ The Flatte parameters are chosen in 
order to reproduce the qualitative features of the data. We find that with Mz c = 3.92 GeV and z = 12 GeV~5 the 
main features of the data can be reproduced. We postpone a more refined extraction of the parameters, including the 
uncertainties, until the tttt final state interaction is included properly. 

As expected, the tttt spectrum of Fig. [5] (b) is in disagreement with the data, however, we expect that this discrepancy 
can be overcome, once we include the 7T7t final state interaction. 

In Fig. [H we also included the production channel for V(4260) — > h c irir. Similar to the discussions of Zb — > h^ir, 
this channel is ideal for disentangling the molecular nature of the intermediate Z c (3900): The power counting analysis 
in Ref. [13] shows that the triangle transition Zb — > hb(mP)ir is not suppressed compared to Zb — > T(nS)Tr, although 
the decay is via a P-wave. This explains why the branching ratios for Z\, — > hb(mP)ir are compatible with those 
for Zb — > T(nS)n. A similar phenomenon occurs here. In addition, higher loop contributions are suppressed in 
Zb — > hb(mP)ir, while they are not suppressed for Zb —t Tn [17]. The analogous pattern is expected for the Z c 
decays. However, as we will see, the box diagrams contribute negligibly to the transition Z c — >• h c ir such that we are 
still able to predict the shape of the h c n invariant mass distributions. Therefore, a measurement of Z c — ¥ h c n and its 
comparison with Z c — > J/ipn will provide a further non-trivial confirmation of the molecular nature of Z c (3900). 

In Fig. [3l we present our prediction for the invariant mass spectra of h c ir and 7rTr in K(4260) — > h c mr. The notation 
of the curves is the same as in Fig. [2] Similar to the J/ipirir channel, a very pronounced peak right at the DD* 
threshold appears in the h c n invariant spectrum. Interestingly, due to the limited phase space in this decay channel, 
its kinematic reflection is significantly shifted and located at higher invariant masses. As mentioned, for the h c ir final 
state the contribution from the box diagram appears to be negligible. Thus, once measured, this channel will allow a 
clean experimental access to the if c (3900) resonance contribution. 

We stress that in order to understand the gross features of the data the diagrams in Fig. Q] can also be replaced 
completely or in parts by tree level diagrams (see, e.g., Refs. 22, 23|). Then, however, the underlying dynamics with 
the interplay of the DD* and DD\ cuts will be lost. We reiterate that it is the DD\ molecule nature of F(4260) that 
provides a natural explanation for the appearance of Z c (3900) in the V(4260) decays, and for other detailed features 
of the spectra. 

In this work we propose that V(4260) is dominantly a _DZ?i(2420) molecule and identify a unique mechanism, namely 
the presence of two-cut condition, which plays an essential role in F(4260) — \ J '/ipmr and h c Trir. We demonstrate that 
without introducing any drastic assumption, the molecular nature of 1"(4260) as a bound state of DDi can naturally 
explain the observation of an enhancement around 3.876 GeV in the J/^tt invariant mass spectrum. The reflection 
of this peak in the invariant mass spectrum produces a lower mass bump around 3.5 GeV. This signature matches 
the experimental data, thus, provides strong evidence for the molecular nature of V(4260). We also demonstrate that 
for a more detailed description of the data it appears necessary to include explicit -Z c (3900) poles, i.e. a resonance 
structure as an isovector partner of V(3872). 

We stress the following important consequences of this prescription that will be reported in a subsequent paper: 

i) Although the nominal Z)Di(2420) threshold is higher than the mass of y(4260), the threshold of DD*n (then 
decaying to DDirir) is much lower. Dominant decays of F(4260) into DD*ir should be regarded as a natural conse- 
quence and would explain the large deficit in width between the total width and its decays into J /i^ttit and J/ijiKK. 
We also predict an asymmetric spectral shape for the decay of V (4260) into DD*tt. This can be clarified by an energy 
scan around the nominal Y mass (analogous to what happens to the shape of V (4660) when being viewed as ip'fo(980) 
molecule 0). 

ii) A possible reason why the isovector partner Z c (3900) of the A"(3872) was not observed in B decays is that the 
V(3872) might be produced via its small cc component which is absent in Z c . Here, Z c (3900) does show up as a result 
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of the proposed molecular nature of Y(4260), which actually allows the S-w&ve DD* pairs to be copiously produced. 

iii) One might also expect a cusp or possible resonance structure at the D*D* threshold in the decay of ^(4415), 
if this state is assumed to be a molecule driven by the nearby D*Dx threshold. 

iv) The exact mechanism at work here might also be the reason for the appearance of the Z\, and Z' b states in T(5S) 
decays, if we assume that T(55) has a sizeable BB\ component. Note that the relative S'-wave threshold for BB\ is 
only about 120 MeV above the T(5S) mass. 
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